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Abstract 

Nano-carbon materials are investigated intensively. In this paper, the edge-state in nanographene 
materials with zigzag edges is studied theoretically. In particular, while the inter-layer inter- 
actions are considered, we prove that edge states exist at the energy of the Dirac point in 
the doubly stacked nanographene, and in the case of the infinitely-wide lower layer case. This 
property applies both for the A-B and A-C stackings. 
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1. Introduction 



The graphite and single layer graphene materials have been studied intensively, since the electric 
field effect has been found in atomically thin graphene films [1] . These materials can be regarded 
as bulk systems. On the other hand, nanographenes with controlled edge structures have been 
predicted to have localized states along the zigzag edges [2] . The presence of the edge states 
have been observed by experiments of scanning tunneling spectroscopy [3,4]. Thus, the studies 
of the edge states are one of the interesting topic of the field. 

Previously, the magnetic switching effect has been found in the process of insertion and ex- 
traction of molecules in activated carbon fibers [5] . The inserted molecules remain in nanometer 
size pores, and give effective pressure to the nanographite clusters. The Pauli susceptibility 
decreases due to the decrease of the magnetic moment magnitude. We have studied the mag- 
netism of the nanographite using the tight binding model including the interlayer interaction 
between neighboring nanographene layers [6-8] . We have found that open shell nature of each 
layer might explain the experimental observations. Recently, interlayer hopping interaction 
effects in stacked graphite have been investigated theoretically [9]. The electronic structures 
change dramatically including the massless Dirac cone and parabolic dispersions, depending on 
the layer numbers from a single graphene, bilayer graphene, to multi layers. 

In this paper, we will study the edge-state in nanographene materials with zigzag edges 
including the inter-layer interactions, extending the work of a single layer [2]. We will prove 
that edge states exist at the energy of the Dirac point in the doubly stacked nanographene, 
and in the case of the infinitely-wide lower layer case. This property can be shown both for the 
A-B and A-C stackings. 

2. Edge states of one graphene layer 

In this section, we review the idea of the edge state which appears along the zigzag line of a 
graphene [2]. Figure 1 shows a grapnene sheet which has one zigzag line at the top most edge 
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of the figure. The graphene extends infinitely in the left, right, and down directions. We will 
constitute a wavefunction of the edge state at the energy £" = as follows. The lattice sites 
are divided in A and B sublattices. The edge atoms belong to the A sublattice. The edge 
state has the amplitude zero at all the sites of the B sublattice. The wavenumber in the one 
dimensional direction is denoted as k, and a is the lattice constant of the unit cell, which is the 
length between the neighboring edge atoms. The condition that the amplitude becomes zero at 
the atom between (n — l)th and nth edge sites, when the nearest neighbor hopping interaction 
t is considered, is 

^ik{n-l)a ^ ^ikna + ^ ^ Q, (l) 

and this gives the amplitude 

X = [-2cos(A;a/2)]e^'=("-^/2)a_ ^2) 

The similar condition, 

^ikna ^ gifc(n+l)a + ^ ^ Q, (3) 

gives 

y = [-2cos(A;a/2)]e*'=("+^/2)a_ 

The condition that the amplitude becomes zero at the sites surrounded by the amplitudes x, 
y, and z is 

X + y + z^O. (5) 

Therefore, we obtain the amplitude at the third zigzag line 

z = [-2cos(A;a/2)]'e^'="". (6) 

Iterating this process, we can prove that the amplitude at the nth unit cell along the mth is 
proportional to [— 2cos(/ca/2)]'"~^. The absolute value is smaller than unity if the wavenumber 
satisfies 27r/3<A;a<7r/2. Therefore, the edge state at the energy E — extends into the bulk 



of the graphene with exponentially damping amplitude with respect to the distance from the 
edge. The edge state locates near the boundary of the first Brillouin zone in the wavenumber 
space. 

3. Edge states in bilayer graphene with parallel zigzag 
lines 

Figure 2(a) shows the A-B stacked bilayer graphene with zigzag edges at the top of the figure. 
The system extends infinitely in the other directions. The upper layer is shown by the solid 
lines, and the lower layer (dotted line) is shift downward by the bond length b. The number 
n ravels the nth unit cell shown by the dashed line. The other label m indicates each zigzag 
line. For the convention, the upper layer begins with m — 0, and the lower layer starts with 
m = 1. At the circles, two carbon atoms overlap completely, and weak hopping interaction ti 
is assigned here. The same m indicates the zigzag lines in the upper and lower layers. 

We will constitute wavefunctions of the edge state. It is assumed that the edge state has 
amphtudes at the sites of the A sublattice only. The amplitude becomes zero in the B sublattice. 
Even if the interactions ti are present at circles, the alternation of the A and B sites remains 
in the whole system. So, the bipartite nature remains, too. When we look at the condition of 
the edge state at the energy £■ = of the upper layer, we find the formula, 

On.me-^'^"/' + an+i.^e^'^"/' + an,m+l = 0, (7) 
for even m, and 

for odd m, where a^^^ is the amplitude at the nth unit cell of the mth zigzag line. Assuming 
Gnft = A, we obtain 

an,m = ^[-2cos(A;a/2)]"^. (9) 
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The similar condition of the zero energy state gives, 

&n,me"^'="/' + 6„+i,^e*^"/' + bn,m+l + rian,m = (10) 

for odd m, and 

for even m, for the amphtude bn,m of the lower layer. Here, ri — ti/tis the ratio of the interlayer 
hopping integral to the intralayer hoppings. As we know the form of a„^^, we can solve the 
recurrence formula of the number series to obtain 

bn,m = -(m - l)nA[-2cos{ka/2)]"'-^ + B[-2cos{ka/2)]"'-^ (12) 

for m > 2, where bi^m — B is assumed. Therefore, we have found that the edge state persists in 
the A-B stack case. The magnitudes of A and B will be nearly equal A B, so the correction 
by the interlayer interaction is of the order ri with respect to the intralayer term. 

Figure 2(b) shows the A-C stacking case, where the lower layer is moved into the right-down 
direction with a bond length b. The notations are similar to those of Fig. 2(a). The difference is 
the fact that the both zigzag edge lines begin with the same index m — 1. In the A-C stacking 
case, it is found that the amplitude of the edge state can be constructed in the lower layer, 
first. The condition of the zero energy state gives, 

bn,me-""'^' + &n+i,me^''^/' + = (13) 

for even m, and 

&n,me-''="/' + 6n+l,me''"/' + &n+l,m+l = (14) 

for odd m, for the amplitude bn,m- This sequence is solved as 

Vm = B[-2cos{ka/2)r-\ (15) 

where bn,i — B. The effects of the interlayer interactions give the next relation for the upper 
layer, 

On.me"'*^"/^ + a„+i,„e^'=«/2 + a„,^+i + ri6„,^ = (16) 



for odd m, and 

On.me"'*^"/^ + a„+i,„e^*^"/2 + a„+i,^+i + ribn+i,m = (17) 

for even m. This is solved with the help of the previous solution bn,m to give 

a„,^ = A[-2cos(ka/2)]'^-^ - (m - l)riS[-2cos(A;a/2)]"*-2 (18) 

for m > 2, where a„,i = A. The roles of the upper and lower layers seem to be exchanged from 
that of the A-B stacking case. The amplitude of the upper layer has correction terms of the 
order ri owing to the presence of the lower layer. 



4. Edge states in bilayer graphene with infinitely-wide 
lower layer 

In this section, we will consider the bilayer graphene, where the lower layer has infinite spatial 
extent. So, the nanographene with one zigzag edge is placed upon one infinite graphene. We 
consider both of the A-B and A-C stackings. They are shown in Figs. 3(a) and (b), respectively. 
As the convention of the choice of the starting point of the initial amplitudes, A and B, is 
different, we would hke to formulate again for the both cases. 

Figure 3(a) shows the A-B stacking case. As the upper layer A-site {n,m = 1) interacts 
directly with the lower layer site as denoted by the circles, the starting points are indexed 
differently from those of the previous section. Here, the zigzag lines of the upper and lower 
layers in the down direction are indexed by m. And, the zigzag lines of the lower layer in the 
up direction are indexed by I with I — m — 1. In the upper layer, the amplitude is 

an,m = A[-2cos(A;a/2)]"^-^ (19) 

with Qn,! — A. In the lower layer, the condition of the zero energy state is the relations, 

i'n.me-^'^"/' + 6n+l,me^''"/' + &n,m+l + nan,m = (20) 
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for even m, 

&n,me"'*^"/' + &n+l,me''="'" + &n+l,m+l + TlO^+i,^ = (21) 

for odd m, 

for even i and I >2, and 

Cn,me-''"/' + C„+i,^e^*^«/' + Cn+l,m+l = (23) 

for odd 1. We obtain, 

bn,m = -mriA[-2cos(A;a/2)]"*-^ + S[-2cos(A;a/2)]"^-^ (24) 

and 

c^i = S[-2cos(A;a/2)]^-^ (25) 

for I > 2, where bn,i — Cn,i — B. 

Figure 3(b) is the A-C stacking case. In the lower layer, the amplitude of the edge state is 
calculated as 

bn,m = B[-2cos{ka/2)Y' (26) 

and 

c„,, = E[-2cos(A;a/2)]\ (27) 

where 6„,o = Cn,o = B. In the upper layer, the effect of the interaction appears. The amplitude 
can be derived by the relation 

On.me-''^"/' + a^+i.^e^'^"/' + a„,^+i + r^hn,m = (28) 

for odd m and m > 1, and 

On.me"'*^"/^ + an+i.^e^'^"/^ + a„+i,„+i + ri6„+i,„ = (29) 



for even m. This is solved with the help of above solution bn,m to give 

a„_^ = A[-2cos(A;a/2)]"*-^ - mriB[-2cos{ka/2)]'^ (30) 

for m > 2, where a„,i = A. Therefore, we have constructed the wavefunction of the zero energy 
state for the infinitely- wide lower layer cases. 

5. Summary 

We have studied the edge-state in nanographene materials with zigzag edges including the 
inter-layer hopping interactions. We have shown that edge states are present at the energy of 
the Dirac point in the doubly stacked nanographene, and in the case of the infinitely-wide lower 
layer case. This property has been found both for the A-B and A-C stackings. 
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Figure Captions 



Fig. 1. The edge states along the zigzag edge of a graphene. See the text for the notations. 

Fig. 2. Bilayer graphene with parallel zigzag lines, (a) A-B stacking and (b) A-C stacking are 
shown. See the text for the notations. 

Fig. 3. Bilayer graphene with the infinitely-wide lower layer, (a) A-B stacking and (b) A-C 
stacking are shown. See the text for the notations. 
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Fig. 1. The edge states along the zigzag edge of a graphene. See the text for the notations. 
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Fig. 2. Bilayer graphene with parallel zigzag lines, (a) A-B stacking and (b) A-C stacking are 
shown. See the text for the notations. 
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Fig. 3. Bilayer graphene with the infinitely-wide lower layer, (a) A-B stacking and (b) A-C 
stacking are shown. See the text for the notations. 
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